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We propose a scheme to dynamically synthesize a space-periodic effective magnetic field for neutral atoms 
by time-periodic magnetic field pulses. When atomic spin adiabatically follows the direction of the effective 
magnetic field, an adiabatic scalar potential together with a geometric vector potential emerges for the atomic 
center-of-mass motion, due to the Berry phase effect. While atoms hop between honeycomb lattice sites formed 
by the minima of the adiabatic potential, complex Peierls phase factors in the hopping coefficients are induced by 
the vector potential, which facilitate a topological Chern insulator. With further tuning of external parameters, 
both a topological phase transition and topological flat bands can be achieved, highlighting realistic prospects 
for studying strongly correlated phenomena in this system. Our Letter presents an alternative pathway towards 
creating and manipulating topological states of ultracold atoms by magnetic fields. 
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Gauge fields lie at the center of our modern understand¬ 
ing of physics in many systems, including those of high en¬ 
ergy and condensed matter, as well as of ultracold atoms. 
Within the gauge-field paradigm, interactions between par¬ 
ticles, which enable rich quantum phases of a many-body 
system, are mediated through gauge fields. For instance, 
solid-state materials with charged quasiparticles in magnetic 
fields or with spin-orbit coupling (SOC) show a rich variety 
of quantum Hall effects and exotic topological superconduc¬ 
tivity [1, 2]. The interplay between gauge fields and lattice 
systems is also of great interest (for a pedagogical review, 
see [3]). The spectrum of a charged particle in a square lat¬ 
tice exposed to a strong uniform magnetic field shows a frac¬ 
tal structure, widely known as the Hofstadter butterfly [4]. In 
another seminal work, Haldane shows that the quantum Hall 
effect without Landau levels can be realized when a periodi¬ 
cally staggered magnetic field is applied to charged particles 
in a honeycomb lattice [5]. 

Ultracold atoms in lattice systems are considered to be pow¬ 
erful simulators for studying gauge-field physics [6-23]. Both 
the Hofstadter and the Haldane models with cold atoms were 
theoretically proposed [11-14] and experimentally demon¬ 
strated [15-21] by making use of novel forms of light-atom 
interactions [7], such as laser-assisted tunneling [15-18], the 
shaking-optical-lattice technique [21], and SOC within a syn¬ 
thetic dimension [19, 20]. In addition to the optical lat¬ 
tice formed from a space-periodic ac-Stark shift by inter¬ 
fering laser beams, proposals for the generation of a mag¬ 
netic lattice with a space-periodic Zeeman shift have been put 
forth [24-32] (and some have been realized [25-28]) using 
current-carrying wires [24], microfabricated wires or perma¬ 
nent magnetic structures on atomic chips [25-30], supercon¬ 
ducting vortex lattice shields [31], and phase imprinting by 
gradient magnetic pulses [32]. In contrast to optical lattices, 
magnetic lattices are free from atomic spontaneous emissions 


that are always accompanied by heating and loss of atoms. 
Additionally, they have the potential to reach shorter lattice 
constants [30, 31] (of a few tens of nanometers as proposed in 
Ref. [31]), leading to improved energy scales and less strin¬ 
gent temperature requirements; the lattice constants can even 
be continuously tuned [32]. These advantageous features en¬ 
hance the performance of atomic quantum gases as powerful 
quantum simulators. 

While the simulation of gauge-field physics and the ma¬ 
nipulation of topological states in optical lattices have shown 
fruitful results [6-23], it remains to show whether this is also 
the case for magnetic lattices. This Letter provides an affirma¬ 
tive first answer to this question. 

This Letter presents a scheme for synthesizing a time- 
independent effective Hamiltonian with nontrivial band topol¬ 
ogy for atomic gases with internal spin degrees of freedom, 
based on the phase imprinting technique [33, 34]. A two- 
dimensional (2D) magnetic lattice with triangular geometry 
emerges in the effective Hamiltonian. In the limit when an 
atom is confined in the lowest Zeeman level, an adiabatic 
scalar potential and a geometric vector potential are simul¬ 
taneously generated for the center-of-mass motion [6, 22, 23]. 
The adiabatic scalar potential surface can form a honeycomb 
lattice, while the associated geometric vector potential pro¬ 
vides complex phases for next-nearest-neighbor (NNN) hop¬ 
ping coefficients in realizing the Haldane model [5, 35]. With 
the flexibility and tunability of magnetic fields, our scheme 
can be extended to produce a set of effective Hamiltonians 
whose lowest energy bands undertake a topological phase 
transition from a topological (Chern) insulator to a trivial one. 
Moreover, models possessing topological quasiflat bands are 
realized near the phase-transition point. 

The protocol .—We consider a pancake-shaped quasi-2D ul¬ 
tracold atomic gas of spin F confined in the x-y plane (at 
z = 0). In the presence of a bias magnetic field BqCz, the 
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single-particle Hamiltonian is given by 

Hq- — + hojoFz, ( 1 ) 

2m 

where p = iPx,Py) is the 2D kinetic momentum, m is the 
atomic mass, h is the reduced Planck constant, is the 
third component of the atomic spin vector (in unit of h) 
F = (Fx, Fy, F^) and wq = EpPeBoITl is the Larmor frequency 
at Bq, where gr is the Lande g factor for the spin-F hyperfine 
state manifold and pB is the Bohr magneton. 

A short gradient-magnetic-field pulse B'ycy of duration 6t' 
imprints a space-dependent phase factor [32-34, 36-38] onto 
the wave function as exp^-iksoyF^,), where k^o — Sf gfpgB' jti 
is the SOC strength [36, 37] with B' the magnetic gradient. 
After a free evolution time 6t, a second magnetic field pulse 
in the opposite direction imprints an opposite phase. The two 
pulses combined together enact a unitary transformation 

eik^yFy F.e~‘'‘«yF ^ F, cos {k,oy) - Fx sin {k,„y ), (2) 

which rotates the magnetic field Bo(0,0,1) to a space- 
periodic form Bq{- sin(ksoy), 0, cos(ksoy)). Similarly, an op¬ 
posite uniform-field pulse pair +Bye^. with a pulse area 
dt'gfpBBylh = 7T inverts the magnetic field Bq (0,0,1) to 
Bo (0,0,-1) as - -F^. More generally, a 

gradient magnetic field pulse along an arbitrary direction 
eg - (cos 0, sin 9, 0) in the x-y plane imprints a phase fac¬ 
tor cxpi-iksofgFg), where rg - r ■ eg - xcosO H- ysint? and 
Fg - F ■ eg - Fx cos 9 + Fy sin 0 are, respectively, the coordi¬ 
nate vector r = {x, y) and the spin vector F projected to the eg 
direction. Following a period of free evolution and a second 
pulse from an opposite gradient field, an expression analogous 
to Eq. (2) generates a magnetic field with spatial periodicity 
along the eg direction. 

In our scheme, discussed below, repeated pulse pairs are 
concatenated. A complete cycle of the evolution period con¬ 
tains three gradient pulse pairs along directions separated by 
an angle of 120°, together with a +7r pulse parr along the y 
direction as shown in Fig. 1(a)-1(b). The total evolution over 
one complete cycle (of period T = 46t) is then given by 

UiT,0) = 

^ih„re.Fe. ^-ihorg.Fg. (3) 

>3,2,1 

with 9j - —nib + 2717 / 3 . According to the Floquet theo¬ 
rem [40, 41], a time-independent effective Hamiltonian can 
be defined according to U {T,0) = cxp{—iH^sTlti). To the 
lowest order of T, we find [42] 

1/3 15 , 

~ 2mF ~ SFPB^eS ■ F, (4) 

where Fj^ = (Fx,Fy) is the 2D spin operator, cjj^q — 1ik^^/2m 
is the SOC frequency, and Beff is an effective magnetic field 



FIG. 1. A schematic illustration for synthesizing a magnetic lattice, 
(a) An atomic cloud is exposed to a uniform bias magnetic field Bge. 
(pointing out of the page), and subjected in sequence to four pairs 
of opposite magnetic pulses as shown in (b). Pulse pairs 1, 2, and 
3 are gradient magnetic fields lying in the x-y plane along directions 
separated by 120°. Pulse pair 4 executes +7r spin rotations along the 
y direction. The color gradient represents the corresponding mag¬ 
netic field strength, (b) The time-periodic pulse sequence with four 
pairs of pulses forms one complete evolution period, (c) A three- 
dimensional view of the effective magnetic field Beu [Eq. (5)], which 
forms a Skyrmion lattice [39] with its Wigner-Seitz unit cell shown 
by the hexagon [of edge length 47r/(3ks(,)]. The arrows are colored 
by the magnitude of the third component of Beff. 

whose three components are given by 

Bes,x = Y^,sin(^ksorg^)sin9j, 

Beft,v = sin [ksoFgj) cos 0j, (5) 

BsS,z = ^ [“1 + {ksoFg^) ■ 

The first two terms in Eq. (4) arise from the unitary transfor¬ 
mations by gradient pulse pairs applied to the momentum op¬ 
erator [32, 36]. The third term describes a magnetic (Zeeman) 
lattice that couples the atomic spin to the effective magnetic 
field Beff, as shown in Fig. 1(c). 

Geometric potentials and energy spectrum. —The above 
protocol for the generation of a triangular magnetic lattice is 
general, and can be applied to atoms with arbitrary spins. For 
concreteness, we choose a specific atomic species, fermionic 
®Fi, with electron spin J - Ijl, nuclear spin 7=1, and we 
consider the total hyperfine spin F — I - J — 1/2 ground-state 
manifold. The Fande g factor can be evaluated according to 
the Breit-Rabi formula [43] tobeg/r -1/3. The spin op¬ 
erator reduces to F = cr/2, where cr is the vector of Pauli 
matrices. To be more specific, in all numerical calculations, 
we assume a set of fixed parameters unless otherwise noted. 
They are Bq = 20 mG, B'5t' = 2 G cm ' ms [44], which corre¬ 
spond to kso = (1.7yum)“' and Wo = 32.3a;so = (27r)x9.3kHz 
for the F - 112 manifold of ®Fi. With these parameters, the 
lattice term in Eq. (4) dominates during time evolution. 

We denote the space-dependent eigenstates of the magnetic 
lattice by [yi, 2 (r)), which satisfy 

gFPB^eS ■ J lTi. 2 (r)> = +eo(r) |^i, 2 (r)>, (6) 
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FIG. 2. Mapping the effective Hamiltonian [Eq. (4)] to the Haldane 
model, (a) The density plot of the adiabatic potential eo (in blue) 
and its associated vector potential A shown by field of arrows. A 
darker color denotes a smaller eo> the minima of which form a hon¬ 
eycomb lattice with two sites per unit cell denoted, respectively, by 
red and green filled circles. The dashed lines denote NNN hopping 
paths along the directions of positive Peierls phases, (b) The flux 
density n^, in units of The hexagons in (a) and (b) denote 

the primitive unit cell [the same as in Fig. 1(c)], over which the net 
flux is unity, (c) Energy spectrum. Solid lines represent the lowest 
two energy bands along lines with high symmetry in the first Bril- 
louin zone (inserted hexagon) for the effective Hamiltonian Eq. (4), 
with (Jo = 32.3^50. Dashed lines represent the fitted band structure 
using the Haldane model results. The edge length of the inserted 
hexagon is VS. (d) A logarithmic plot of the Berry curvature for 
the lowest band fli. The integration of Di in the first Brillouin zone 
gives its Chern number, Ci = 1. 


where eo(r) = is the adiabatic potential for 

atomic center-of-mass motion in the lower-energy eigenstate 
\X\)- For an atom adiabatically moving in this space-periodic 
Zeeman level, a vector potential A(r) emerges [42, 45], 

3 

A = y^i) H-;t^soCrilcr± in)- (7) 

lo 

with (Tj^ = (cTj., <jy). Associated with the vector potential is 
the flux density = (V x K)^l2nh, which shares the same 
spatial periodicity as Befr and can be considered in general as 
a type of flux lattice [22]. 

The adiabatic potential eo, vector potential A, and the flux 
density for our magnetic lattice are shown in Figs. 2(a) 
and 2(b) [46]. As shown in Fig. 2(a), the local minima of eo 
are located at the corners of the unit cell, forming a honey¬ 
comb lattice. When an atom hops between these honeycomb 
sites, the vector potential contributes a complex Peierls phase 
factor exp^/ J A ■ dl/h^ to the hopping coefficient [4, 5, 35], 
with the integration evaluated along the corresponding hop¬ 
ping path. As A vanishes along the edges of the hexagon, the 
nearest-neighboring (NN) phase factor is a trivial unity. While 
along the NNN hopping paths [dashed lines in Fig. 2(a)], the 
accumulated phases are always nonzero. Thus the adiabatic 
scalar potential together with the geometric vector potential 


realizes the Haldane model in the tight-binding limit. As a 
caveat, our flux pattern shown in Fig. 2(b) is not the same as 
that suggested by Haldane [5], where the staggered flux den¬ 
sity gives a vanishing net flux over a unit cell. The flux dis¬ 
tribution shown in Fig. 2(b) is non-negative everywhere, and 
the net flux over one unit cell is unity rather than zero; this 
can be checked by integrating the following over a unit cell: 

£fuc ^ ^ Beff/|Beff| [22]. 

Thus, the nontrivial winding pattern of Bgff shown in Fig. 1(c) 
leads to a quantized net flux = 1. A nonzero net flux gen¬ 
erally leads to larger Peierls phases (of order unity). It also 
facilitates simulation of charged particles in strong magnetic 
field with nondispersive Landau levels [22, 47]. 

To quantitatively confirm that our model indeed maps onto 
the Haldane model, we numerically study the spectrum and 
Berry curvature [48] of the effective Hamiltonian Eq. (4) us¬ 
ing the plane-wave expansion method [42, 49]. The typical 
band structure and the Berry curvature for the lowest band are 
shown, respectively, in Figs. 2(c) and 2(d). A band gap opens 
at the corners of the first Brillouin zone {±K points), where 
the Berry curvature is at a maximum. Both the eigenenergies 
and the Berry curvatures are even functions of quasimomen¬ 
tum, so the spectrum at K' — -K is not shown. The Chern 
numbers [50] for the lowest two bands are Ci _2 = ±1, respec¬ 
tively. The spectrum and the Berry curvature thus resemble 
the ones from the Haldane model. To further validate this cor¬ 
respondence, we adopt the method used in Ref. [51] to get 
the NN hopping constant H and the complex NNN hopping 
constant |f 2 k"^ of the Haldane model from the calculated band 
structure. We find H = 0.053^Wso and |f 2 | = 0.0037^Wso with 
4> - 0.40. Using these three parameters together with an over¬ 
all energy shift, the tight-binding band structure of the Hal¬ 
dane model is plotted with dashed lines in Fig. 2(c). 

Topological phase transition and quasifiat bands .—Our 
protocol allows for the easy tuning of two parameters: the 
SOC strength k^o - df gfpsB'Ih and the bias magnetic field 
Bq. Both can be tuned continuously, and can be turned on 
adiabatically to reach the ground state for our model system 
Eq. (4) [8, 21] (see [42] for details). Once the ground state 
is achieved, we can apply an additional weak optical gradient 
field (which commutes with all the pulse manipulation oper¬ 
ations) in the x-y plane to drive Bloch oscillations and then 
measure the perpendicular center-of-mass drift to extract the 
topological properties for the lowest energy band [17, 21, 52]. 
With unequal durations between subsequent pulse pairs, or al¬ 
lowing for specific kso and Bq values for different subperiods, 
several variants of the effective Hamiltonian can be synthe¬ 
sized. A topological phase transition for the lowest energy 
band can be achieved by a simple tuning of the bias magnetic 
field. Eor this to occur, we set the field strength to be Bq for the 
first three subperiods and switch to aB^ for the fourth subpe¬ 
riod; our protocol then leads to a change for the z component 
of the effective magnetic field in Eq. (5) as 


Bq 


-a + 


Xi,cos(^sore,) 


( 8 ) 
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FIG. 3. Illustration of the structural and topological phase transition 
with the tuning of the z component of B^jf in Eq. (8). (a) Top panel: 
the density plots of the adiabatic potentials for a = 0 (left), a = 
0.720 (middle), and a = 1 (right) with the minima (dark colors) 
forming a simple triangular lattice, a decorated triangular lattice, and 
a honeycomb lattice, respectively. Bottom panel: The Chern number 
for the lowest energy band as a function of a. A topological phase 
transition occurs at the critical point of O' = 0.720. (b) The band 
gap between the lowest two bands and the band width for the lowest 
band, as a function of a. The band gap closes at o- = 0.720 (marked 
by the dotted vertical line). The inset shows the dependence of the 
band gap-over-width ratio on a. It peaks at o- = 0.749 with a value 
of 16. (c) The lowest three bands at tt = 0.749. 


The a - 1 case corresponds to the original proposal with topo¬ 
logical bands, while the cr = 0 case describes a system of triv¬ 
ial energy bands with zero Chern numbers. By continuously 
tuning a from 0 to 1, a topological phase transition with band 
touching and reopening takes place, as summarized in Fig. 3. 

Figure 3(a) presents the changing Chern number and, 
hence, the band topology, for the lowest band with increas¬ 
ing a. The lattice geometry of the adiabatic potential is found 
to undergo a structural transformation first from a simple tri¬ 
angular lattice to a decorated triangular lattice [53], and then 
to a honeycomb lattice. The corresponding tight-binding de¬ 
scriptions for s orbitals involve 1, 3, and 2 bands, respectively, 
for the three cases. As a increases, the band originating from 
hopping between s orbitals located at unit cell centers crosses 
the two bands from s orbitals located at the corners. Their 
corresponding Chern numbers change after band touching and 
reopening. Figure 3(b) shows the behavior of the gap between 
the lowest two bands as well as the band width for the lowest 
one. Gap closing occurs at the F point when a - 0.720, and 
the gap-over-width ratio is found to be quite large over a lim¬ 
ited range after a gap opening with a peak value as large as 16 
when a - 0.749, as shown in the inset of Fig. 3(b). The energy 


spectrum at cr = 0.749 is shown in Fig. 3(c). The lowest band 
is a Landau-level-like topological quasiflat one [54]. Such a 
nondispersive topological band also persists beyond the adia¬ 
batic limit [42]. It is a promising candidate for simulating the 
fractional quantum Hall effect when suitable interactions are 
included [55-58]. It is perhaps worth pointing out, also, that a 
flat band can emerge as the second excited band in a Kagome 
lattice [53], or as the first excited band in a Lieb lattice [59]. 
The properties of the localized states in the flat band of a Lieb 
lattice have been investigated in a recent experiment [59]. 

We focus in this Letter on discussing single-particle physics 
of a fermionic spin-1/2 system, though our magnetic lattice 
generation protocol can be equally applied to a higher-spin 
atom, be it a boson or fermion. When local momentum- 
independent (s wave) interaction is taken into account, it 
can be simply added to the effective Hamiltonian because 
it commutes with all the pulse manipulation operations (see 
also [36, 37]). The topological phase is expected to be sta¬ 
ble for weak interactions due to the presence of a gap. How¬ 
ever, stronger interaction can drive the system to new phases, 
in which the physics may be dominated by the interplay be¬ 
tween the correlation and band topology. A detailed study of 
the interaction effects for this system deserves further efforts. 

In conclusion, we propose an experimentally feasible pro¬ 
tocol to realize a synthetic magnetic held with real magnetic 
field pulses. The synthetic magnetic field forms a lattice 
with nontrivial band topology, and under certain limits can 
be mapped to the Haldane model. The high tunability of our 
scheme makes it possible to design a topological phase transi¬ 
tion as well as quasiflat energy bands with nontrivial topology, 
which could push the effective model into the strongly corre¬ 
lated regime. 
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This supplementary material provides more details for the various points mentioned in the main article. 


Derivation of the effective Hamiltonian 


As explained in the main text, a gradient magnetic field pulse along eg - (cos 0, sin 0, 0) in the x-y plane prints a phase factor 


exp (-iksofgFg) onto the single particle wave function. Under the action of an opposite pulse pair, the atomic momentum p = fik 
and spin transform according to. 


^ _ I^^^Pg Sm0, 


(SI) 


where Fg = F^ sin 0 - Fy cos 0 is the spin operator perpendicular to eg. Thus the evolution including the pulse pair gives rise to 



(S2) 


(S3) 


A +7t pulse pair along the y direction flips the direction of the magnetic field seen by atom without inducing coupling to 
momentum. The evolution operator from such a pulse pair is given by 


Uy,A^t,0) = 


(S4) 


with the corresponding Hamiltonian 



(S5) 


A complete evolution period contains three gradient pulse pairs separated by a mutual angle of 120° in the x-y plane, together 


with a +7T pulse pair along the y direction, the corresponding evolution operator for one period T - A6t is given by 

U (T - 4(5f, 0) = 


= g-iHcsT/n 

with 0j - We use to denote the Hamiltonian for each subperiods as 


(S6) 


(S7) 


Hi = Hg^, H2 = Hg,, Hi = Hgy, H^ = Hy,, 


To lowest order of the Trotter’s expansion for T, the effective Hamiltonian //eft is the average of 



(S8) 
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with = 4 for the current case. Expanding out the terms in this Hamiltonian explicitly, we arrive at the effective Hamiltonian 
as shown in the main text after recombining them. 

The higher order corrections of the effective Hamiltonian can also be derived through the Magnus expansion method, which 
to first order of T gives [41] 


“ At Xi 2nN^- ^ 


(S9) 


m<n=2 


From Eq. (S6) we see that, our protocol can also be viewed as a periodically driven 4-step sequence, where the dynamics 
of the driven system can be understood in terms of the effective Hamiltonian and the associated micromotion [8]. We now 
provide the details for how our system fits into this formulism. The Hamiltonian for a general N-step sequence is described by 
+ Vm, with the constraint 2m=i - 0. Thus we find 




which is just the time-averaged effective Hamiltonian Eq. (S8), and correspondingly 

1 ^ 


(SIO) 


(Sll) 


m=l 


The evolution operator for such a driven system can be partitioned in the following as introduced in Ref. [8] 

U (t, 0) = (S12) 

where the effective Hamiltonian and the initial-kick operator K(Q) of the N-step sequence to first order of T are given by [8] 


Heft = = H^°'> Yj 




m<n=2 


(S13) 




with Cm,n - j + m -n. For a full evolution cycle (t - T), the evolution operator is found to be 

U(T 0) ^ e^'K(T)^-iTH,sin^iK{0) ^ ^-iTH,g/n 

where the Hamiltonian H^g to first order of T is given by [using K{T) - /r(0)] 

Hgg = ^ -H - /[Z:(0),//™]. 


(S14) 


(S15) 


Inserting Eqs. (SIO), (Sll) and (S13) into the above equation, we recover Eq. (S9) as expected. 

This shows that, the micromotion of the N-step periodically driven system gives the first order correction to the effective 
Hamiltonian. This calls for the use of the time-averaged effective Hamiltonian Eq. (S8) to approximate the evolution in the short 
T limit. This approximation will be validated numerically in the later sections. 


Solving the effective Hamiltonian with the plane wave expansion method 


For spin-1/2 with F = (t/ 2, the effective Hamiltonian takes the following form 


1 


2m 


Hgg = — p - ■ a-. 


16 


1 


(S16) 


with (Tj^ = cTy) after omitting a constant energy shift 15Swso/128. As H^g in Eq. (S16) is space-periodic, its eigenstates can 
be labeled with quasimomentum q = i^qx, qr) in the first Brillouin zone as good quantum numbers. Expanded in the plane wave 
basis [49], the eigenstates take the form 




'(/bi+mb2)-r 




(S17) 
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where ii is the band index, bi _2 = ^so (+cos - sin are the two reciprocal unit vectors, l,m are integers taking values 
0, +1, +2, etc. The expansion coefficients are to be determined. The plane wave expansion Eq. (S17) together with the 

eigenvalue equation 


(r) = E„ (q) ^„q (r), 


(SI 8) 


determines the energy spectrum, as well as the Bloch wave functions as eigenfunctions. 

In our numerical calculations, we use a cutoff = 15 as the maximal value for |/| and \m\. And we have checked that, for 
A > 10 cases, the band structures shown in the main text are independent of the choice of N. Thus in the plane wave basis, the 
Hamiltonian is expressed as a [2 (2N + 1)^] x [2 (2N + 1)^] matrix for each quasimomentum q within the first Brillouin zone. 
The eigenvalues and corresponding eigenvectors give respectively the energy spectrum and Bloch wave functions. 

Once the Bloch wave functions are obtained, we can calculate their Berry curvatures Q„ (q) and (first) Chern numbers C„ for 
each energy bands according to [48] 


^11 (tj) ~ ^ 

c„ = 7^ I d-qn„ (q), 

27r Jbz 


(S19) 


where M„q(r) = (r |M„(q)) = e '‘^''^iA„q(r) is the cell-periodic part of the Bloch function. 


Validity of the effective Hamiltonian 

We use the effective Hamiltonian to describe the evolution of ultracold atoms under the magnetic field pulse sequence of the 
proposed protocol. The use of Trotter expansion limits its validity to short evolution periods. In this subsection we use two 
complimentary methods to check for this approximation. 



0 0.2 0.4 0.6 0.8 1 

time (ms) 


FIG. S1. The evolution of a Gaussian wave packet located at the center of a harmonic trap, with oscillator frequencies = co,. = (2;r) x 30 Hz 
is investigated. The system is initially prepared in the spin up state. The figure shows the time dependent population imbalance (cr.) calculated 
from the effective Hamiltonian (solid line), and the actual pulse sequence (open circles), with evolution periods T = 40ps, 20ps, and lOps, 
from top to down respectively. The other parameters are cuo = 32.3tUso = Q-n) x 9.3 kHz, a = 

The first method replies on evolving an eigenstate, whereby the initial state i/^„q (r, f = 0) is prepared in an eigenstate of the 
effective Hamiltonian with a given band index n and quasimomentum q. The actual magnetic field pulse sequence is then 
used to evolve it to a later time f > 0 where the wave function is denoted by (^„q (r, t). The overlap between these two states 
(t) |iAnq (t = 0)) should be identically equal to unity at all times for the ideal case when the effective Hamiltonian exactly 
represents the evolution by the pulse sequence. 

We arbitrarily choose a quasimomentum point q = (0.3, -Q.2)kso in the first Brillouin zone, and the eigenstate for the lowest 
band n = 1, for the case too - 32.3tUso = (27r) x 9.3 kHz and a - 1, is calculated according to the method described in the last 
section. The absolute values of the overlap after evolving for 1 ms for T - 40 ps, 20 ps and lOps, are respectively calculated to 
be 0.947, 0.984 and 0.995, indicating improved level of approximations with shorter periods. 
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The second method relies on wave packet evolution. For this case, an external trapping potential Vaa^{x,y) - \m{Lo\x^ + 
is included into the Hamiltonian, and the initial state is prepared as a wave packet in the harmonic trap. Since the 
trap potential commutes with all the pulses, it is simply added into the effective Hamiltonian. The initial wave packet is then 
evolved respectively through the effective Hamiltonian and through the actually pulse sequence. We can then calculate some 
physical observables, e.g., population imbalance, by the corresponding time-dependent states, and compare their respective 
results. Figure SI shows an example of this comparison. We can see that, they coincide with each other better for shorter 
evolution periods. 

These studies indicate that for sufficiently short evolution periods, the effective Hamiltonian faithfully describes the evolution 
under the actual magnetic field pulse sequence. For the cases we discuss, T = 10/rs < 1 /wo = 17/zs is short enough to validate 
the effective Hamiltonian approximation. 


Derivation of the emergent gauge fields 


The evolution of a single-particle state is governed by the effective Hamiltonian according to 

iS-|'T(r,t)) = //efFl'I'(r,t)), (S20) 

ot 

with //eff given by Eq. (S16) for the spin 1/2 case. Expanded in the adiabatic basis, with corresponding center-of-mass wave 
functions labeled as if/j (r, t), the wave function takes the form 

|'T(r,t))= 2^,(r,t)^,(r)), (S21) 

> 1,2 

which when acted upon by the spin-dependent shift leads to 

(p - IT' (r, t)) = 2 - Ab- - Ay) Wi) , (S22) 


with 


Ay = in {xi\ , 

3 

Ay = <Xi\ \xj) ■ 


(S23) 


In the adiabatic limit wq » Wso, if the initial state is prepared in the dressed state \x\), then the probability amplitude for the 
particle to be in the orthogonal state \x 2 ) remains zero at all time. Thus by projecting Eq. (S20) to the dressed state \xi) and 
taking ^2 - 0, we get a closed equation for i/fi 


^ d 
ot 


—(p - A)^ + ea + W 
Zm 


iAi> 


where 

3 

A = All + All = in(xi\ Vxi) + —nk^o{Xi\<^i. iTi)> 

fo 

is the geometric vector potential that couples to the center-of-mass motion, and 

W - -— (Ai2 ■ A 21 + Ai2 ■ A 21 -I- Ai2 ■ A 21 -t- Ai2 ■ A 21 ), 
2m '■ ’ 


(S24) 


(S25) 


(S26) 


is the geometric scalar potential. 

The resulting scalar potential consists of two terms. One is the adiabatic potential eo, whose characteristic energy scale is htOQ. 
The other is the geometric scalar potential W, which scales as natso- Thus W makes a negligible contribution to the total scalar 
potential when compared to eo, in the wo » Wso limit. In Pig. S2, geometric potentials for a = 0, 0.72, and 1 at wo = 32.3wso 
are shown. Compared to the corresponding adiabatic potentials shown respectively in Pig. 3(a) in the main text, with the 
corresponding potential (minima, maxima) in units of /iWso being (-12.0, -4.0), (-9.2, -6.7) and (-10.0, -7.0) respectively, the 
geometric potentials only give small corrections to the final total scalar potentials. The qualitative understanding of the Haldane 
model as well as the topological phase transition enabled by our model is not affected by including these corrections. 

If we take Ay = 0 in Eq. (S23), then the above results reduce to those reviewed by Dalibard et al. in Ref. [6], where momentum 
is not coupled to spin components in the original Hamiltonian, and the corresponding geometric potentials are A = iH (yJ Vyi) 
and W = Ai 2 ■ A 2 i/ 2 m = C^il V;r 2 ) lV2m. 
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FIG. S2. The geometric scalar potentials W (in units of ftw^o) in Eq. (S26) at = 32.3Wso for (a) or = 0, (b) cr = 0.72, and (c) a = 1. 


Energy spectrum beyond the limit of adiabatic approximation 


Although the discussion in the main text on the generation and understanding of the topological energy bands is based on 
models under the adiabatic as well as tight-binding approximation, the nontrivial topology of the energy bands is found to exist 
beyond these two approximations [22, 47]. As an example, we take a - 0.749, and reduce the original wq value by a factor of 
five to 0)0 = = (27r) x 1.9 kHz. The term of a Zeeman level is used instead of the adiabatic potential when the adiabatic 

approximation is not satisfied. At ojq - 6A6a>so with a = 0.749, the (minima,maxima) of the two Zeeman levels of the magnetic 
lattice are respectively given by (-1.84,-1.34) and (1.34,1.84) in units of ^Wso- The separation between the corresponding two 
levels is 2.68?zWso, and the lattice depth of the lower Zeeman level is a mere O.S^Wso. As the energy scale for the spatial uniform 
term in Eq. (S16) is of order Ha>so, and SOC in this term can flip the spin, the lattice term does not dominate during the evolution; 
the adiabatic approximation fails for this case. 



FIG. S3. The energy spectrum for our periodically driven model system at o = 0.749 and oiq = 6.46tJso = (In) x 1.9 kHz. (a) From the 
effective Hamiltonian Eq. (S16). (b-d) The quasienergies from solving the evolution operator in Eq. (S6) for (b) T - 50 /js; (c) T = lOOps; 
and (d)r = 200ps. 


The lowest six energy bands of the effective Hamiltonian for this case are shown in Fig. S3(a). Their Chern numbers are 
all found to be €„ - 1, which invalidates directly a tight-binding description, because the sum of Chern numbers for a set of 
tight-binding bands should be equal to zero. The band width for the lowest energy band is 0.0049/iWso, and the lowest band gap 
is Q.22tKx)f,Q, which gives a gap-over-width ratio of ~ 45. Thus these Chern numbers for the lowest few bands resemble the ones 
for a charged particle in a magnetic held in a weak lattice background [49]. 

Next we discuss the validity of the effective Hamiltonian for this case. In previous subsection and in Fig. SI, we show the 
validity of the effective Hamiltonian for wq = 32.3wso = {2.n) x 9.3 kHz, with respectively evolution periods T - lOps, 20ps, 
and 40 ps. In the present case of loq - {2ji) x 1.9 kHz, i.e., reduced by a factor of five, the suitable evolution period is expected to 
be increased by a similar factor. We can repeat the previous method by scanning state evolution period to illustrate the validity 
of adiabatic approximation, or alternatively, we can calculate the Floquet spectrum of quasienergies by evaluating the evolution 
operator Eq. (S6) directly for different evolution periods, with the results as shown in Fig. S3(b-d). We see that the spectrum of 
the quasienergies strongly resembles the one from the effective Hamiltonian at least for T - lOOps. The topological properties 
of the lowest energy band is expected to be preserved for T - 200 ps, because no signature of band closing is found. For even 
larger evolution period, e.g. T - 400ps, the quasienergies for the lower energy bands are found to be smeared out due to 
the wrapping of the higher Floquet sectors with quasienergies at multiples of 2ntilT. Thus we conclude the suitable evolution 
period required for the validity of the effective Hamiltonian approximation is of the order of ~ lOOps for the present choice of 
0)0 = (27r)x 1.9 kHz. 










































11 


Adiabatic preparation of the ground state 

As is verified in the previous sections, one can use the effective Hamiltonian Eq. (S16) to describe the state evolution under 
the pulse sequence for short evolution periods. In this section, we describe how can the ground state of the effective Hamiltonian 
Eq. (S16) be reached through an adiabatic loading approach. Eor a non-interacting many-body fermionic system, the ability of 
preparing its ground state is essentially equivalent to the statement that all the (single-particle) eigenstates in the lowest energy 
band of in Eq. (SI 6 ) can be reached appropriately. 

As shown in the main text, one complete cycle of the evolution period consists of four subperiods. To prepare the ground 
state, we will keep the fourth subperiod intact, while adiabatically ramp up (a) the strength of the gradient pulses in Na complete 
cycles and then (b) the bias magnetic field for the first three subperiods in Nt, complete cycles. 

A quantitative analysis is provided as follows. We first introduce an auxiliary Hamiltonian 

3/jLn / \ a y 

- - -+ Py^y) ~ ■ cr, (S27) 

2m lorn y ''8 8 

where the three components of the M vector are respectively - - Yjj ^in {k^org^ sin (Oj^, My - Yjj sin (ksorgj^ cos ^0^^, and 

= Yjj cos The effective Hamiltonian before the aforementioned ramping process is ’H (k = Q,y - 0); the effective 

Hamiltonian at the end of the ramping stage (a) is TY (/c = 1,7 = 0); and the effective Hamiltonian after ramping stage (b) is 

(k - 1 , 7 = 1), which recovers the effective Hamiltonian //eft in Eq. (S16). 

As a proof-of-principal illustration, suppose the initial state for the ramping stage (a) is prepared as a plain wave in the spin 
up state 


I *Aini,a) 


= exp (ip ■ r/K) 


vOy 


(S28) 


which is the eigenstate of "K (a- = 0 ,7 = 0) in the lower branch (we always consider positive a). The target state for the ramping 
stage (a) is the corresponding eigenstate of "K (/c = 1,7 = 0 ): 


<r I lAtar.a) = eXp (/p ■ T/H) 


3uJsoiPx + iPy)IKo 




(S29) 


Here the spin state is not normalized. We denote the wavefunction during stage (a) under the pulse sequence as |0p (f))- The 
fidelity between this state and the target state is then given by 


We notice that, the initial fidelity for stage (a) is 


Ep (t — 0) — l^^tar.a | ^ini.a^j 


which is close to unity for small momentums and large crcuo/wso- For instance, for the momentum point p = ( 1 , \ )1ik^o in the 
lower energy branch, with wo = 32.3wso = (27r) x 9.3 kHz and a - 1, we get Ep (f = 0) = 99.57%. After linearly ramping 
up K from zero to one in Na - 40 cycles with period T - 10/rs, the final fidelity is numerically found to be Ep (f - NaT) - 
1 - 4.3 X 10“^. This shows that the low energy eigenstates of TY (a- = 1,7 = 0) can be reached reliably. 

We then consider the ramping process for stage (b). We choose the initial state for this stage as the lowest eigenstates of 
Ef(A' = 1 , 7 = 10“^). The vanishingly small lattice term 10 ® x 4^woM ■ cr introduced in our numerical simulation accounts for 
the effect of the Bragg reflections at the edges of the first Brillouin zone, which turns the good quantum number from momentum 
(denoted as p) to quasimomentum (denoted as q) restricted in the first Brillouin zone [49]. The fidelity between the target state 
l^iq^ [i.e., eigenstate of Eq. (S16) with the form of Eq. (S17) with band index n = 1] and the time-dependent state |0q (f)) 
evolved by pulse sequence with increased magnetic lattice strength (which is proportional to 7 ) is defined as: 

Eq(f) = |(lAlq I <^q(0)| ■ 


' ~ I^^Atar.a | 0p ■ 


(S30) 


iatiLOQ + ^J9htOsoP^/2m + (cr/zwo)' 


iahoJo + -^9/jWsoP^/2ot + (ahcoof] + 9/jWsoP^/2ot 


(S31) 


(S32) 
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We suggest ramping up the strength of the magnetic lattice as a tangent function (rather than a linear one), which corresponds to 
a time-dependent j function with the form 


r(0 = 


tan \riTitj {2NbT)\ 
tan (rjnll) 


(S33) 


where the parameter 77 is taken as 0.975. The typical results for the time-dependent fidelity functions for different quasimomen- 
tums, together with the y{t) function, are shown in Fig. S4. As can be seen in this figure, the fidelities after ramping up the 
magnetic lattice approach unity for all tested quasimomentum points. We thus conclude that the eigenstates in the lowest energy 
band of the effective Hamiltonian //eff in Eq. (SI 6 ) can be prepared appropriately through the adiabatic approach we propose. 



FIG. S4. The time-dependent fidelities for selected quasimomentum points in the first Brillouin zone. Here in the legend, P stands for the 
result for the quasimomentum point q = (0.3,-0.2)kso. The strength of the magnetic lattice is ramped up as a tangent function (the inset) in 
Nb = 1000 cycles. The parameters for the system are the following: tuo = 32.3 £Js(, = (In) x 9.3 kHz, a = 1 and T = lOps. 














